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SUMMARY 


The role of Ames Research Center in conducting basic aerodynamics 
research through computer simulations is described. The computer 
facilities, including supercomputers and peripheral equipment, 
represent the state of the art. The methodology of computational 
fluid dynamics is explained briefly. Fundamental studies of tur- 
bulence and transition are being pursued to understand these 
phenomena and to develop models that can be used in the solution 
of the Reynolds-averaged Navier-Stokes equations. Four applica- 
tions of computer simulations for aerodynamics problems are 
described: subsonic flow around a fuselage at high angle of at- 
tack, subsonic flow through a turbine stator-rotor stage, tran- 
sonic flow around a flexible swept wing, and transonic flow 
around a wing-body configuration that includes an inlet and a 
tail . 


INTRODUCTION 


Since the introduction of jet aircraft to regular commercial 
service 30 years ago, the world's fleet has grown to over 7,000 
airplanes. Although the United States' market share has eroded 
over the years, most of these aircraft were manufactured in the 
United States. New aircraft are needed now to replace aging 
airframes, to comply with noise restrictions, and to improve 
operating efficiencies. Manufacturers plan to build over 8,000 
new aircraft during the next 15 years. One of NASA's objectives 
is to provide the aircraft industry with the basic technology re- 
quired to design new aircraft that will maintain the country's 
preeminence in the world marketplace. 


COMPUTER FACILITIES 


Ames Research Center has two computer facilities available 
for the simulation of aircraft aerodynamics. The Central Comput- 
ing Facility serves the general research needs of the Center with 
a Cray Y-MP/832 supercomputer. The Numerical Aerodynamic 
Simulation (NAS) facility provides a computational capability to 
the whole country, including researchers in industry, univer- 
sities and government agencies. The high-speed processors in the 
NAS facility are currently a Cray-2 and a Cray Y-MP/832. 
Thirty-five scientific workstations are distributed among users 
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at the Center; other workstations are distributed around the 
country. The complete system includes mass storage, support 
processing, and long-haul communications subsystems connected by 
a high-speed data network. The long-term goal is to replace the 
high-speed processors periodically with the fastest available su- 
percomputers . 


COMPUTATIONAL FLUID DYNAMICS 


The Navier-Stokes equations which describe the behavior of 
viscous flows have been known for over 150 years. Computational 
fluid dynamics is the relatively new field in which these partial 
differential equations are solved using high-speed computers. 
The flow field, including any solid bodies or boundaries, is en- 
cased within a grid of mesh points. Difference methods are then 
used to solve the flow properties for the specified boundary and 
initial conditions. 


TURBULENCE 


Nearly all practical aerodynamics problems involve turbulent 
flow near surfaces and in wake regions. The combination of the 
fine grid required to resolve the turbulence and the processing 
speed required to perform the calculations rules out exact solu- 
tions of the Navier-Stokes equations for complex configurations. 
In order to obtain solutions, the turbulence terms are time- 
averaged and approximated with empirical models. 

The most commonly used turbulence model is the Baldwin-Lomax 
algebraic eddy viscosity model (Ref. 1), which works well for 
flows with little or no separation. Improved models are neces- 
sary for flows with large separated regions and steep pressure 
gradients. Computer simulation of turbulence is being used to 
study the detailed structure for the purpose of developing ne / 
turbulence models. Spalart (Ref. 2) recently performed direct 
simulations of a turbulent boundary layer on a flat plate for 
Reynolds numbers, Re 0 , based on momentum thickness, up to 1410. 
The detailed quality of the calculations may be observed in the 
vorticity contours shown in Fig. 1. 


APPLICATIONS 


Current capabilities for computer simulation of aerodynamics 
problems are demonstrated by the following applications made by 
members of the Ames Applied Computational Fluids Branch. The 
equations solved are usually the thin-layer, Reynolds-averaged 
Navier-Stokes equations, and the Baldwin-Lomax model is used for 
turbulent flow. 
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Fig. 1 Vorticity contours on flat plate, Reg = 1410 (Ref. 2) 


^ingraf t , fore b e fl y 

Subsonic flow around a slender fuselage at a high angle of 
attack was calculated by Schiff, Cummings, Sorenson, and 
Rizk (Ref. 3) to study the vortical flow on the leeward side, 
which affects the maneuverability of the aircraft. A closeup of 
the two-block grid is shown in Fig. 2; the interface is located 
at the beginning of the wing leading-edge extension (LEX) . The 
aft block has a larger number of circumferential points to define 
the LEX. Nearly a quarter of a million grid points were used. 



Fig. 2 Two-block grid around Fig. 3 Surface flow pattern on 
aircraft forebody aircraft forebody 

(Ref. 3) (Ref. 3) 

The surface flow pattern is shown in Fig. 3 for M«, = 0.2, 
a = 30°, and Reynolds number based on mean aerodynamic chord, 
Re c = 11.5 x 10 . The primary and secondary separation lines on 
the forebody and the secondary crossflow separation line on the 
LEX are in good agreement with flight-test data. 

T u r b ine s tage 

Subsonic flow through a turbine stator-rotor stage was cal- 
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culated by Rai (Ref. 4) to understand the flow processes in order 
to design more efficient turbomachinery. A patched-grid system 
was used with the grid about the rotor moving relative to the 
grid about the stator. This is shown in Fig. 4 for a midspan 
location. The complete geometry included the hub and outer 
casing with some tip clearance. Equal number of rotor and stator 
blades were used to reduce the calculations to one rotor and one 
stator blade. The total number of grid points was 200,000. 


ZONE 3 



Fig. 4 Patched grid around 
stator and rotor at 
midspan (Ref. 4) 



Fig. 5 Instantaneous pressure 
contours at midspan 
(Ref. 4) 


Calculations were performed for an inlet Mach number, 
M^ = 0.07, and Re = 10 5 per inch. Instantaneous pressure con- 

tours at midspan are shown in Fig. 5. 

Fl ex ible w ing 

Transonic flow around a flexible swept wing was calculate 
by Guruswamy (Refs. 5,6) to demonstrate the coupling between the 
vortical flow and the wing deformation. The unsteady Navier- 
Stokes equations were solved simultaneously with the structural 
equations of motion. A new grid was generated at each time step 
to conform with the deformed wing. The initial and deformed 
grids at the 50% semispan airfoil section are shown in Fig. 6. 

Calculations were performed for ramp motion of a rectangular 
wing, using a NACA 0015 airfoil section with an aspect ratio of 4 
and a sweep angle of 30°. The effects of flexibility on un- 
steady lift coefficients at four span locations are shown in Fig. 
7 for Moo = 0.5, Reynolds number based on chord, Re c = 2 x 10 , 
and pitching rate A = 0.1. The higher lift for the flexible 
wing is attributed to the increase in angle of attack caused by 
the flexibility. 
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Fig. 6 Initial and deformed 

grid at 50% semispan of 
flexible wing (Ref. 5) 


Fig. 7 Unsteady lift coeffi- 
cient for rigid and 
flexible wing (Ref. 6) 


Complete aircraft 

Transonic flow around a wing-body configuration, including 
an inlet and a tail, was calculated by Flores and Chaderjian 
(Ref. 7). The purpose was to demonstrate the ability to calculate 
transonic viscous flow over a complete aircraft. The flow field 
was divided into 27 zones, determined by geometrical or flow con- 
ditions, for a total of half a million grid points. The surface 
grid covering the forward portion of the fuselage, including the 
inlet and wing, is shown in Fig. 8. 

Computations for = 0.9, a = 6°, and Reynolds number 
based on wing root chord, Rec = 4.5 x 10 , required 25 hours for 
5,000 iterations, using one processor of a Cray X-MP/48. Pres- 
sure coefficients on the vertical tail at three-quarters height 
are shown in Fig. 9. Comparison with experimental results 
(Ref. 8) shows good agreement. Calculations for the wing and 
fuselage without the tail surfaces (Ref. 9) also show good agree- 
ment (Fig . 10) . 
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8 Grid for fuselage 
forebody (Ref. 7) 


Fig. 9 Pressure coefficient on 
vertical tail (Ref. 7) 



10 Pressure coefficient on wing and upper fuselage center- 
line (Ref. 9) 


CONCLUDING REMARKS 


Many challenges still remain in computational fluid 
dynamics. Steady-state solutions are sought using time-dependent 
methods that require thousands of steps for convergence. The ef- 
ficiency of numerical methods must be improved to make the best 
use of available computers. New turbulence models are needed; 
the most commonly used model — the Baldwin-Lomax model — is over 
ten years old and is inadequate for separated flows. Single- 
processor supercomputers are approaching their inherent perfor- 
mance limits, and multiprocessors offer a means for increasing 
performance. Algorithms and codes must be constructed to take 
advantage of the supercomputer architecture. Finally, aerody- 
namics must be combined with other aspects of aircraft design. 
Structure and materials, control and guidance, and propulsion 
systems are equally important, and coupling between disciplines 
will be required in future designs. 
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